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The extensive application of surfactants motivates comprehensive and predictive theoretical studies 
that improve our understanding of the behaviour of these complex systems. In this study, an 
expression for the elastic free-energy density of a wormlike micellar chain is derived taking into 
account interactions between its constituent molecules. The resulting expression incorporates the 
sum of a quadratic term in the curvature and a quadratic term in the torsion of the centerline 
of wormlike micelle and thus resembles free-energy density functions for polymer chains and DNA 
available in the literature. The derived model is applied on a wormlike micelle in the shape of 
a circular arc, open or closed. A detailed application of the derived model on wormlike micelles 
of toroidal shape, along with employing necessary statistical-thermodynamical concepts of self- 
assembly, is performed, and the results are found to be consistent with the ones available in the 
literature. Steps towards obtaining the material parameters through experiments are suggested and 
discussed. 


I. INTRODUCTION 

A surfactant molecule consists of two main parts: a 
hydrophobic tail and a polar hydrophilic head-group. 
When present in solution at sufficiently high concen¬ 
trations, surfactant molecules self-assemble into vari¬ 
ous supramolecular structures that shield the hydropho¬ 
bic tails from contact with the ambient solution (Is- 
raelachvilpJ21) These structures include spherical mi¬ 
celles, short cylindrical micelles, long cylindrical micelles 
called wormlike micelles, bilayers, and closed bilayers 
or vesicles (Cates and CandarP and Discher et a/.®). 
A collection of spherical micelles may undergo uniaxial 
growth and coalesce into short rod-like cylindrical mi¬ 
celles (Cui et alW and Groswasser et The ends of 

a cylindrical micelle are capped by hemispheres (Porte 
et al^. Adding more molecules at certain temperatures 
and osmotic pressures leads to the formation of worm¬ 
like micelles (May and Ben-ShauP and Dam et 
The energy required to create two end-caps from a 
very long cylindrical micelle is called the scission energy 
(Oelschlaeger et If this energy is sufficiently large 

and the volume fraction of the surfactant molecules is 
sufficiently low, the semiflexible micelles may fuse to min¬ 
imize the number of end caps (Cates and CandaiP). 

In recent years, considerable attention have been paid 
to investigating the rheological properties of surfactant- 
based micellar solutions (Boek et al.^ Lerouge and 
Berret,^^ and Spenley et Some examples include 

soap solutions including micellar globules and adhesives 
formulated with block copolymers (Won et al^^. Due 
to their novel properties, such solutions have found ap¬ 
plications as heat-transfer fluids, hard-surface cleaners, 
liquid dish-washing detergents, drag-reducing agents in 
pipelines, and fracking fluids (Cecile,^ Yang,l^ and Ko- 
mura and SafraiP^. 

Wormlike micelles may be regarded as mesopolymers 
since they impart elasticity (Padding and BoeMiS and 
BerretP^ and are subject to breaking and re-shaping 


(Boek et al.^^ Zhou et al.^^ and Germann et 
They may also become entangled or form branched struc¬ 
tures (Andreev and Victorov^Sl). As a result, a solution 
in which they are suspended m ay ha ve viscoelastic be¬ 
haviour (Acharya and Kunieda,ESl2ll Raghavan et alW^ 
and Yan^i^. For this reason, many rheologists focus 
on understanding the effective non-Newtonian behav¬ 
ior of micellar solutions (Kuperkar et and Rad- 
ulescu et alW^. That behavior is generally influenced by 
micro-structural changes that occur during flow (Kha- 
tory et The present paper focuses on the ener¬ 

getics of an individual closed or open wormlike micelle, 
leaving aside questions related to dissipative interactions 
in flowing solutions. 

Wormlike micelles are often characterized by their per¬ 
sistence length. This quantity can be obtained from 
the bending free-energy (Nettesheim and Wagner,!^ and 
Schubert et al^^. The free-energy density therefore has 
a significant role in studying these materials. May et 
considered the extent to which the bending elasticity of 
wormlike micelles influences their tendency to join and 
form branched structures. In their work, the free-energy 
density of a wormlike micelle is comprised of the chain 
conformational free-energy, the end-cap energies, and the 
hydrocarbon-water interfacial energy. Their findings in¬ 
dicate that the energy change associated with the forma¬ 
tion of a junction between one micellar end-cap and the 
cylindrical body of another micelle is small relative to the 
free energy of an end cap. 

The first proposed conformational free-energy density 
function for a worml ike micelle was motived by the 
C anham-HelfricfP^*^ elastic free-energy density for a 
lipid vesicle, 

ijj = ^o + ]^kc{H - Hof+ kcK, ( 1 ) 

where 70 is the surface tension, H and K are the mean 
and Gaussian curvatures of the surface, kc and kc are the 
splay and saddle-splay moduli, and Ho is the spontaneous 
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mean curvature, namely the mean curvature of the natu¬ 
ral, local shape of the bilayer that has been attributed to 
the difference between the volumes of the head and tails 
(Tang and Carteil^. 

Lauw et used self-consistent field-theory to esti¬ 
mate the free-energy of a closed toroidal wormlike micelle 
comprised of nonionic surfactant molecules. This relied 
on specializing the Helfrich expression Q to a toroidal 
geometry. Proceeding similarly, BergstronP^^^ investi¬ 
gated the effect of the splay modulus on the size and 
shape of toroidal micelles and their stability. 

An alternative approach to express the free en¬ 
ergy of a wormlike micelle is based on the well- 
established molecular-statistical perspective for model¬ 
ing self-assembled aggregates. In such an approach, the 
free energy per amphiphile is considered as the sum of 
various contributions including hydrophobic effect, chain 
conformational entropy, electrostatic effects (obtained 
from linearized Poisson-Boltzmann equation), and head- 
group repulsion (Puvvada and Blankschtein,!^ Nagara- 
jan and Ruckenstein,ISl and May et According 

to Bergstrom,!^ the aforementioned contributions can 
be accurately described by applying Helfrich’s curvature 
elasticity theory 0 . 

Wormlike micelles have diameters on the order of 3-5 
nanometers. The contour length of a wormlike micelle 
need not be fixed and may be sensitive to environmen¬ 
tal conditions. The length of a wormlike micelle may 
be extremely large relative to its cross-sectional diame¬ 
ter, reaching 3600 nano meter s (Cecil^^^ or as much as 
one millimeter (Paddin^i^Illl) j-atio of the length 

of a polymer chain to its thickness can be up to 500 
or more (Roiter and MinkcP^. The same ratio for a 
typical human hair ranges between 10^ and 10^. The 
length to diameter ratio of a wormlike micelle exceeds 
the aforementioned ratios for polymer chains or human 
hair. Considering the geometric analogy between worm¬ 
like micelles and polymer chains, and the application of 
rod-like models for pol ymer chains (Kratky and Porod,!^ 
and Bugl and FujitcP^, an approach closer to those used 
to model polymers seems to be worthy of consideration. 

As previously mentioned, the notion of free-energy 
density is central to the study of wormlike micelles. For 
instance, with an expression for the free-energy density 
of such a micelle, equilibrium configurations can be ex¬ 
plored. Here, a free-energy density function is derived 
for wormlike micelles. To do so, an approach based on 
accounting for the interactions between the constituent 
molecules of the micelle is adopted. This approach is 
largely based on ideas developed by Keller and Mer¬ 
chant who showed that the surface energy of a sub¬ 
stance generally includes a bending term and, thus, that 
the surface may exhibit elastic resistance to bending even 
in the absence of stretching. This contribution to the free 
energy is obtained in terms of the molecular density and 
interaction potential. In a recent application of the work 
of Keller and Merchant ,I^Seguin and FriecP^ derived the 
Canham-Helfrich free-energy density for a lipid vesicle. 


In so doing, they modeled the lipid molecules comprising 
the vesicle by one-dimensional rigid rods. For simplicity, 
they neglected molecular tilt relative to the orientation 
of the bilayer and interactions between the bilayer and 
the solution. 

To determine the free-energy density of a wormlike mi¬ 
celle at a position x, we account for the interactions be¬ 
tween all surfactant molecules within a cutoff distance d 
from the molecules at x. For simplicity, we model sur¬ 
factant molecules as one-dimensional rigid rods and as¬ 
sume that these rods are perpendicular to the centerline 
of the wormlike micelle. Additionally, interactions be¬ 
tween the surfactant molecules comprising the wormlike 
micelle and the surrounding solution are neglected. Our 
derivation relies on a Taylor series expansion with re¬ 
spect to a dimensionless parameter S := dji <C 1, where 
I is the smallest radius of curvature that the centerline 
of a wormlike micelle is capable of exhibiting. The free- 
energy of an open wormlike micelle results on integrating 
this density over the centerline of the wormlike micelle 
and adding the end-cap energies. 

This article is organized as follows: In Section |TTj 
modeling assumptions and geometrical considerations are 
briefly presented. The section is followed by introduc¬ 
ing the distribution function of the molecules along the 
body of the wormlike micelle. Section m includes the 
calculation of the free-energy density along the tubular 
body of the wormlike micelle, followed by obtaining the 
free-energy of closed and open wormlike micelles. Sec¬ 
tion |IV| is concerned with the application of the model 
to a wormlike micelle in the shape of a circular arc, open 
or closed. A detailed application of the model on worm¬ 
like micelles of toroidal geometry along with statistical- 
thermodynamical concepts and suggested experimental 
procedure towards obtaining the material parameters are 
provided in Section]^ Finally, the key findings and im¬ 
plications of this study are summarized and discussed 
in Section lYg Details of the various derivations are pre¬ 
sented in the Appendix. 

II. MODELING ASSUMPTIONS 

According to the assumptions below, the body of a 
wormlike micelle is envisioned as a tubular domain with 
centerline C, of length L, and constant cross-sectional 
radius a: 

1. The surfactant molecules that comprise the micelle 
have identical physiochemical properties. 

2. Each surfactant molecule can be modeled as a one¬ 
dimensional rigid rod of length a. 

3. The surfactant molecules at any point of C are ar¬ 
ranged perpendicular to C with a uniform angular 
distribution. 

4. The length L of the wormlike micelle far exceeds its 
radius a, namely the length a of a single molecule. 
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5. The minimum radius of curvature that a wormlike 
micelle can support at each point is denoted by i. 


Moreover, the curvature k and torsion r of C are given 

by 


As a consequence of assumptions and the head- 
groups of surfactant molecules lie on a tubular surface of 
constant circular cross-section. Assumptionis based on 
the observation that the hydrophobic tails of surfactant 
molecules are oriented (on average) along the normal of 
the tubular interface (May et and Shikata et . 
Assumption allows the wormlike micelle to be identi¬ 
fied with its centerline C. Assumption refers to the 
fact that the tubular body of the micelle may adopt var¬ 
ious shapes; however, being made up of molecules of a 
finite size, it cannot support arbitrarily large curvatures. 
Moreover, the tubular body does not overlap (or even 
contact) itself. On this basis, there is an upper bound 
for the curvature (or, equivalently, a lower bound for the 
radius of curvature) which ensures that contact/overlap 
does not occur. We denote such radius of curvature by 
i. 


K = It'l = lx" 


and T = ±|b'| = 


x' • (x" X x'") 


■ ( 6 ) 


Let s belong to the interval (0, L), so that x(s) is in¬ 
terior to C. Consider a molecule at the position x(s) 
with orientation 9 measured counterclockwise from the 
line determined by the unit normal n(s). The director of 
this molecule is given by the unit vector d(s, 9). Accord¬ 
ing to assumption]^ the director d(s, 0 ) of the molecule 
at x(s), can be expressed as a linear combination 


d(s, 9) = {cos 9) n(s) - 1 - (sin0) b(s), (7) 

of unit normal n(s) and unit binormal b(s) (Figure]^. 
Similarly, the director e{t,r]) of the molecule at position 
x(t) can be expressed as a linear combination of n{t) and 
h{t) by (Figure]^ 

e{t, T]) = (cos 77 ) n(t) -f (sin 77 ) h{t). ( 8 ) 


A. Preliminary geometrical considerations 


B. Molecular distribution function 


Let C denote the centerline of the tubular domain that 
represents the wormlike micelle and consider an arclength 
parametrization 

C = {x : X = x(s),0 < s < L}. (2) 


Using a prime to denote differentiation with respect to 
the arclength s, it follows that 



/dx dx\i/2 
V ds ds/ 


and, also, that 


x' • x" = 0 , and |x' x x"| = |x"|. 


( 3 ) 

( 4 ) 


The unit tangent t, unit normal n, and unit binormal b 
of the Frenet frame {t, n, b} of C are given in terms of 
the arclength parametrization x of C by (Figure]^ 


The way that surfactant molecules are distributed 
along the tubular body of the wormlike micelle, has 
an effect in describing how such molecules interact. 
The molecular distribution function describes how such 
molecules are distributed along C. 

The distribution of the molecules at a generic point s 
in the open interval (0,L) along C is denoted by / > 0. 
It follows that the integral 

f{s) d9 ds (9) 

represents the total number of surfactant molecules com¬ 
prising the tubular body of the wormlike micelle. In view 
of assumption in Section [IT] the angular distribution of 
surfactant molecules alone C is uniform. Consequently, 
([^ yields 

27r [ /(s)ds. ( 10 ) 

Jo 




FIG. 1: Two positions x(s) and x(t) on the centerline of a wormlike 
micelle with the intermolecular vector r, and the Frenet frame at 
those positions. 


III. FREE-ENERGY OF A WORMLIKE 
MICELLE 

Wormlike micelles are either closed or open. An open 
wormlike micelle possesses two end caps. Branched 
wormlike micelles, however, possess more than two end 
caps along with junctions which may contribute addi¬ 
tional free energy (Jain and Bate^^ and Dam et 
Here, we first derive the free-energy density on the body 
of the wormlike micelle. The net free-energy of a closed 
wormlike micelle is simply found by integrating the de¬ 
rived free-energy density over the body of the micelle. For 
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FIG. 2: Geometry of a part of the tubular body of a wormlike 
micelle based on simplifying assumptions. The head groups of the 
molecules at x(s) must lie on the boundary of a disk (of radius a) 
in the plane spanned by n(s) and b(s). 


a distance d may interact. On denoting the relative posi¬ 
tion vector between a pair of molecules at x(s) and x(i) 
by r = x(s) — x(i), it follows that 17 satisfies 

|r| > d ^ r2(x(s), x(t), d(s, 9),e{t, 77 )) = 0. (12) 

It is assumed that the cutoff distance d is small relative 
to the minimum radius of curvature t introduced in as¬ 
sumption!^ so that d i or, 

5-.= j« 1. (13) 

Granted that G is frame indifferent, it must obey 

G(x(s),x(t),d(s,6»),e(t,77)) 

- IrP 

= • d(s,6»),r • e(t, 77 ), d(s, 77 ) • 6 ( 1 , 77 )), 

(14) 


an open wormlike micelle, the free energy corresponding 
to the end caps should be added to the net free-energy of 
the body. For simplicity, we restrict attention to worm¬ 
like micelles that do not possess branches and, thus, pos¬ 
sess only two end caps. 


A. Free-energy density of the tubular body of the 
wormlike micelle 


where a factor of 2 has been introduced to simplify sub¬ 
sequent calculations. The quantity G on the right-hand 
side of (14) is a function of four scalar arguments. These 
arguments include the length of the intermolecular vector 
r, and the dot products r • d and r • e between the direc¬ 
tors and that vector and the dot product d • e between 
the directors. Such dot products are related to the angles 
ai and a 2 between the relative position vector r and the 
directors d and e, and also to the angle between the 
two directors d and e, by 


The interaction energy between two rod-like structures 
is based on three ideas, the first of which entails consid¬ 
ering the molecules as one-dimensional rigid rods, the 
second embodies the notion that the rods are made up 
of material points that interact with each other, and the 
third being the principle of material frame-indifference, 
which states that the constitutive relations describing the 
internal interactions between the parts of a system should 
not depend on the external frame of reference used to 
describe them. This principle places restrictions on the 
constitutive equations (Truesdell and Noll^. 

Assume that the interaction between two molecules is 
governed by an energy function depending on the loca¬ 
tion of the two molecules an d the directors indicating the 
orientation of the moleculesl ^^ l ^^ l Consider two molecules 
at positions x(s) and x(t) on C, as depicted schematically 
in Figure From 0 and Q, their directors are of the 
form d(s,^ and e{t,ri), respectively. Let 


f-d = cosai, f-e = cosa 2 , d-e = cosQ; 3 . (15) 


with f being the unit vector corresponding to the inter¬ 
molecular vector r. According to (12), G depends implic¬ 


itly upon d through S, while G does not depend upon d. 
The term d'^ in the denominator of the first argument on 
the right-hand side of (14) is introduced to remove that 


dependence. It follows from (12)-(14) that G satisfies 


s > i G(s^,s(5cosai,s(5cosa2)COSQ;3) = 0. (16) 

The total free-energy Etot of the tubular body of the 
wormlike micelle is given by 

nL pL ^277 p27T 1 

Etot= / / / ^f^(x(s),x(t),d(s, 6 »),e(t, 77 )) 

JO JO JO JO ^ 

f{s)f{t) d6dr]dtds. (17) 


G(x(s),x(t),d(s, 6 »),e(t,? 7 )), (11) 

denote the interaction energy between these two 
molecules. An interaction energy of the form may 

encompass different effects such as steric or electrostatic 
interactions between surfactant molecules. Following 
Keller and Merchant ,1^ assume that the interaction en¬ 
ergy between two molecules separated by a distance 
greater than some fixed cutoff distance d vanishes. To ex¬ 
press this differently, assume that only molecules within 


To reiterate, G(x(s), x(t), d(s, d), e(t, 77 )) is the interac¬ 
tion energy between a molecule at position x(s) with ori¬ 
entation d(s, 0 ) and a molecule at x(t) with orientation 
6 ( 1 , 77 ). Whereas /(s) gives the density of the molecules 
at x(s), f{t) gives the same quantity at x(l). By inte¬ 
grating from 0 to 27r with respect to the variables 9 and rj, 
the interaction energy between all of the molecules of the 
two cross-sections at x(s) and x(l) is taken into account. 
The integrand in ( [T7| ) is multiplied by 1/2 because oth¬ 
erwise it double counts the interaction energy between 
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each pair of molecules. By integrating over t from 0 to 
L, the interaction of all the molecules of the other cross- 
sections with the molecules located in the cross-section 
at the position x(s) is incorporated. The second integra¬ 
tion, over s from 0 to L, ensures that the free energy of 
each cross-section is counted. 

The net free-energy function i?tot of the tubular body 
of the wormlike micelle is related to the free-energy den¬ 
sity Ip by 

Etot = [ V' ds. (18) 

^0 


By @-([ 13 , the free-energy density ip of the wormlike 
micelle at a position x = x(to) on C is, up to an arbitrary 
additive constant, given by 


L p277 p27r 


ip = 


fi 




,r(t) • d{to,0),r{t) ■ e{t,r]), 


lo Jo Jo '• 

d(to,6») • e(t,77))/(<o)/(t) d0dr]dt, (19) 


where r(t) = x(to) — x(t) is the intermolecular vector. 

Upon performing the Taylor expansion of the right- 
hand side of (19) up to two derivatives, the specific steps 


of which appear in Appendix, the hnal form of the free- 
energy density for a wormlike micelle is found to be 


Ip = ipo + kiK^ + A:2T^ 


( 20 ) 


which includes the sum of a quadratic term in curvature 
K and a quadratic term in torsion r of the centerline C 
of the micelle, defined in Q. Notice that the free-energy 
density tp in (20) is independent of the size and shape 


of the wormlike micelle. Thus, according to our theory, 
the free energy of a wormlike micelle in an equilibrated 
system is completely determined by the parameters tpo, 
ki and ^ 2 . 


The term ipo in (201 is insensitive to the shape of the 


wormlike micelle. Since the molecular distribution func¬ 
tion / has implicit dependence upon effects like tempera¬ 
ture, concentration, and electromagnetic fields, these ef¬ 
fects may be encompassed in tpo and in the moduli fci and 
k 2 ■ The role of the term ^po in (1 m I is similar to that of 70 
in Helfrich’s free-energy density (^. However, 70 in Q 
is not an independent parameter, and is determined by 
the condition of constant overall surfactant concentration 
(BergstronP^. 



FIG. 4: Schematic of a wormlike micelle that possesses torsion. 


The flexural (or bending) rigidity ki in (20) represents 
the resistance of the wormlike micelle against deviations 
from a uniform curvature. It further describes the stabil¬ 
ity and stiffness of the wormlike micelle (GradzielskP^l). 
Similarly, the torsional rigidity ^2 denotes the resistance 
of the wormlike micelle against moving out of the osculat¬ 
ing plane of the centerline C at each point. Notice that 
the moduli fci and k 2 depend upon the size and com¬ 
position of surfactant molecules comprising the micelle, 
the distribution of the molecules, temperature, and the 
concentration of the solution (BergstronP^. 

The torsion r of C should not be confused with the 
notion of the twist between different cross-sections of C. 
Since the molecules comprising the wormlike micelle may 
adjust their positions in response to the twist between 
adjacent cross-sections, as shown in Figurej^ the expres¬ 
sion (19) for the free-energy density ip contains no con¬ 
tribution related to twist. In contrast, r measures the 
turnaround of the binormal b of C or, equivalently, the 
tendency of C to move out of a given osculating plane 
(Figure]^. Thus, r vanishes for plane curves and is con¬ 
stant for helical curves. As defined in Q 2 , t is positive 
(negative) for a right-handed (left-handed) helix. 

The expression in (20), which is known in the context 
of the elasticity of bent rods (Landau and Lifshit , re¬ 
sembles free-energy density functions for polymer chains 
(Kratky and Porod,I^Bugl and Fujita,l^and Liu et alW^ 
and those for DNA (Marko and Siggia,l^^l^ and Balaeff 
et However, in contrast to those models, (20) con¬ 

tains evidence of neither intrinsic curvature nor intrin¬ 
sic torsion. For a wormlike micelle consisting of more 
than one type of surfactant molecule, the presence of 
molecules with different head-group or tail conformations 
might lead to intrinsic curvature or intrinsic torsion. In 
the present setting, such effects are ruled out by our first 
modeling assumption. 



FIG. 3: Schematic of a small section of a wormlike micelle. Since 
the molecules can rotate freely, no energetic cost in incurred by 
twisting. 


B. Free energy of a closed wormlike micelle 

For a closed wormlike micelle, the centerline C is closed. 
Integrating the free-energy density function in (20) over 
C yields the total free-energy. 


Etot = J {ipo + kiK^ + fc 2 T^) ds. 


of a closed wormlike micelle with centerline C. 


(21) 
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C. Free energy of an open wormlike micelle 


Consider now an open wormlike micelle with two end 
caps. The net free-energy of such a wormlike micelle is 
obtained by adding the free energy of the tubular body 
and that of the two end caps. The former is obtained by 
integrating the free-energy density (20 1 over C. The latter 
is obtained by 2 times the free energy ifcap corresponding 
to an end cap. The result is given by the sum 



FIG. 6: Toroidal wormlike micelle with constant radius of 
curvature R. 


Etot — J (V'o + kiK^ fc 2 T^) ds -|- 2ii/cap- 


( 22 ) 


B. Example 2: toroidal wormlike micelle 


IV. ILLUSTRATIVE EXAMPLES 

In this section, we investigate the application of our 
model on two examples of wormlike micelles, the first of 
which an open wormlike micelle in the shape of a circular 
arc, and the second, a closed wormlike micelle in the 
shape of a torus. Bearing in mind the relatively high 
difference between the free energy associated with the 
end caps and that associated with the tubular body of 
the micelle, we then explore the condition under which 
the open wormlike micelle in the shape of the circular arc 
tends to close itself to form a toroidal micelle. 


A. Example 1: open wormlike micelle 


Consider a planar open wormlike micelle in the shape 
of a circular arc of length L = R'j, where R and 7 re¬ 
spectively denote the radius of curvature and the central 
angle corresponding to the arc. Assume that the length 
of the arc is large enough, compared to the cross-sectional 
diameter of the wormlike micelle, to ensure assumption]^ 
in Sectionjl^ Since the torsion r of a circular arc vanishes, 
the corresponding net free-energy Etot in ( | 2 ^ simplifies 
to 


Etot — ( ('0O -b ds -|- 2ifcap- 
Jo 


(23) 


Assuming il)o and ki to be uniform along C, while bearing 
in mind that the curvature k is 1/R along the arc and 
ds = AO/k = i?d0, (23) specializes to 


Etot — 


ki \ 
R^J 


2Er. 


(24) 



FIG. 5: Schematic of a planar open wormlike micelle in the 
shape of an arc of a circle. 


Consider a closed wormlike micelle in the shape of a 
torus with the constant major radius R, and the minor ra¬ 
dius a, as indicated in Figure]^ Notice that the length of 
such a micelle, which is L = 2t:R, is assumed to be large 
enough relative to the minor radius a, so that assump¬ 
tion]^ in Section]^ holds. Considering the fact that tor¬ 
sion r for the centerline of a torus vanishes, and k = l/R 
along the centerline C of the torus, while keeping in mind 
that ds = Rd9, the total free-energy Etot in (21) takes 
the form 


k^ 

^^tot = y^ {iPo + j^)Rd0. 


(25) 


Assuming tpo and fci to be uniform along C, (25) results 


Etot = 2nR 


R^J' 


(26) 


Granted that tpo and ki are positive, the free energy 
functions expressed in (24) and (26) tend to -boo when 


R approaches zero or -boo. Also, the second derivative of 
the functions expressed in (24) and (261 with respect to 


R are positive. Thus, both of the functions are convex 
with respect to R, and each possess a single minimum. 
As a result, the corresponding equilibrium state for each 
case is stable. Notice that the minimum of the total 
free-energy Etot in (26) occurs at R* = ^/kij^. Thus, 
our theory predicts that toroidal micelles with constant 
radius of curvature tend to have the major radius R* at 
the equilibrium state. 


C. A comparison between the two examples 

It has been reported that the free energy due to the 
end caps is relatively high (Gelbart and Ben-ShauP^and 
In et al^^. This relatively high free-energy might be 
due to the way that the surfactant molecules are packed 
at the ends to form two semi-spherical caps, which shield 
the water from the hydrocarbon chains (Gelbart and Ben- 
Shauf^. In equilibrium, the free energy of a collection of 
wormlike micelles is minimized by reducing the number of 
end caps (IsraelachvilP). Under certain conditions, this 
may be accompanied by elongation of wormlike micelles 
(Sharma et al^^. 
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Consider the two examples discussed in sections IV A 


and IV B As it transpires, at a specific value jcr of 7 , 
the total free-energy Etot of the open wormlike micelle 
in the shape of the circular arc in (241 becomes greater 


than the total free-energy Etot of the toroidal micelle 
in (26 1 . At equilibrium, the system prefers the state with 
the lower free-energy. Therefore, for 7 > -jcr, the open 
wormlike micelle is expected to display a tendency to lose 
its end caps and form a toroidal micelle; alternatively, 
the micelle prefers to remain open if 7 is smaller than 
7 cr- The angle jcr at which this occurs, can be found by 
equating the right-hand sides of (24) and ( 261 : 


( 27 r - 7 cr)(^ + V’oT?) = 2 Acap- 


(27) 


Hence, 


2R Acap 
7cp = 27r - “P 


ki -I- ijjoR^ 


(28) 


V. APPLICATION OF THE MODEL ON 
TOROIDAL MICELLES AND COMPARISON 
WITH PREVIOUS STUDIES 


In this section, we expand the result of Section IV B 
by employing neces sar y statistical-th ermodyn amical con¬ 
cepts (Israelachvili,lU2l and BergstroirPHHHUl) ^ niake a 

comparison between the results of our model and the pre¬ 
vious studies. 

On thermodynamical grounds, the total free-energy 
Etot of the process of self-assembly of N surfactant 
molecules into a single toroidal micelle can be considered 
as the sum of two terms: the interaction free-energy of 
forming a toroidal micelle out of N surfactant molecules, 
which we denote by Nip-^iic^ and the unfavourable free- 
energy AGtor of self-assembling surfactants, which is the 
positive expression AGtor = —TAStor (BergstronpH) . It 
follows from the set of thermodynamics equilibrium con¬ 
ditions (BergstronP^ 


AGtot — Af'0mic + AGtor — 0, 


(29) 


that ' 0 mic should be a negative quantity in order for the 
process to be thermodynamically feasible. The ideal en¬ 
tropy S'free of mixing N surfactant molecules with 
solvent molecules is 


-S'free = -ksiN In -f Vg In ((> 5 ), (30) 

where (f>s and (phee are the volume fractions of the solvent 
and surfactants respectively. Similarly, the entropy Stoi- 
of the mixture of a single toroidal micelle in Ng solvent 
molecules is 


-Star = -fcs ( In ())Ar -I- Vs In (/)s) , 


(31) 


where is the volume fraction of a single toroidal mi¬ 
celle (BergstronPSl) . Considering (30) and (31), one can 


simply obtain the entropy change AStor corresponding 
to the self-assembly of N surfactant molecules to a single 
toroidal micelle as 


AStor — Stor Sfree — k^ (^ln (jyjy V In . (i^ 2 ) 


In view of (32) and 

AGtor = -TAStor = fcBT( In - V In c/jf ,,,), (33) 

the right-hand side equality in ([ 2 ^ simplifies to 


VV'mic + ksT {hl(pN - N\n c/lfree) = 0 . (34) 


On introducing the net free-energy Atot of the formation 
of a toroidal micelle out of N surfactant molecules as 
(IsraelachvilPI) 


Etot ■= Vf/'mic - V/csTlnf/ffree, 


(34) can be expressed alternatively as 


Ttot + fcsTlnt/fAT = 0. 


(35) 


(36) 


Consequently, 


(/)Ar =exp(-^^). (37) 

ksT 

The total volume fraction i^tot may be expressed 
as the summation of the size distribution function 

(BergstroirPH) 


V^tot 




(38) 


which, according to Bergstrom,*^ can be approximated 
by 



(39) 


In (39), r denotes the ratio of the major radius R of 
the torus, to its minor radius a (see Figure]^. Follow¬ 
ing Bergstrom,!^ the total volume fraction density $tot 
is the integrand of the right-hand side of ( [^ . Hence, 
the total volume fraction density (or the size distribution 
function) <I>tot of a toroidal micelle comprised of N sur¬ 
factant molecules is given in terms of its net free-energy 
Etot{r) by 


^ dN ( Etot{r)' 


(40) 


The quantity N can be obtained by equating the vol¬ 
ume of the tails of N surfactant molecules comprising 
the micelle, and that occupied by the toroidal micelle 
(BergstroirpSl). On representing the volume of the tail 
of a single surfactant molecule by u, the volume of N 
surfactant molecules forming the toroidal micelle is Nv. 
On the other hand, the volume of a torus with the minor 






















FIG. 7: Schematic of the dimensionless net free-energy of a toroidal micelle against the dimensionless radius r of the torus, 
based on our theory; (a) The value of the flexural rigidity was set to ki/a = hksT (BergstronP^, and some arbitrary values 
for ijja = aipo/ksT were selected; (b) The parameter ipo was set to 1, and the value of ki/a was selected between 1-20 ksT. 
The convexity of energy plots indicates the stability of the corresponding equilibrium state. Similar to the previous reports 
(Bergstrorn^, positiveness of the flexural rigidity is a necessary condition for stability. 




r = R/a 


FIG. 8: Schematic of the volume fraction density <I>tot of toroidal micelles in terms of r = R/a, based on our model. Arbitrary 
factors have been used to make each set of the plots in a single frame; (a) The value of the flexural rigidity was set to 
ki/a = 5 ksT (BergstrorrP^ and the parameter ■i/jo = ai/o/ksT was selected between 0.01-1; (b) The value of was set to 1, 
and the flexural rigidity k\/a was selected between 2-10 ksT. 


radius a and major radius R can be obtained by applying 
the second Pappus-Guldinus theorem.l^ Hence, 

Nv = {2'K){'Ka?)R = 2TT^a?r, (41) 


which results 


dN _ 2 -k'^o? 
dr V 


Consequently, (40) becomes 
27r^a^r 


4>t,nt = 


exp 


£^tot(r') \ 

keT )' 


(42) 


(43) 


According to our theory, the net free-energy i?tot of a 
toroidal micelle is expressed in (26). Considering r = 
R/a, (26) takes the form 


Etot{r) = 2-Kar{^o 




fci 


(44) 


or, equivalently, the dimensionless form 


Etotjr) 

ksT 



(45) 


Figure displays the dimensionless free-energy 
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Etot{r)/kBT of a toroidal micelle against the ratio 
r. Granted that 'ijjo and ki remain positive, the net 
free-energy _Etot is convex and possesses a minimum 
point corresponding to the stable equilibrium state. 
Further, the range of the obtained values of the min¬ 
imum free-energy of a toroidal micelle in Figure is 
on the same order of magnitude with that obtained by 
Bergstrom.!^ 

Applying ( [4^ in (431 results the final form of the vol¬ 
ume fraction density ‘f’tot- The schematic of $tot for a 
set of previously reported values of the flexural rigidity 
ki , and some arbitrary values of the parameter ^/>o is de¬ 
picted in Figure As Figure [Sja indicates, size distri¬ 
bution function of toroidal micelles decreases by increas¬ 
ing the value of the flexural rigidity fci, which is con¬ 
sistent with the previous findings (BergstrorrP^. The 
histogram of experimental data can be fitted to the the¬ 
oretical size distribution function $tot in to find the 
real values of the material properties fci and ipo- Recall 
that the size distribution function $tot is a function of 
the total free-energy iftot, which depends on -ipo and fci 
of the toroidal micelle. To simplify yet provide physi¬ 
cally meaningful curve fit, by fixing the value of flexural 
rigidity fci , one can obtain the corresponding best values 
of ipo by fitting the size distribution function <l>tot over 
the histogram of the experimental data. Specifically, one 
can rely on a range of representative values for the flex¬ 
ural rigidity of wormlike micelles reported in previous 
studies. For instance, Jung et alW^ reported the flexu¬ 
ral rigidity in the range of 1-8 fcsT for cetyltrimethy- 
lammonium bromide ^TAB) -based aggregates. More 
specifically, BergstrorrpSl set the flexural rigidity of a sin¬ 
gle toroidal micelle made of ionic surfactants to be 1-5 
ksT. Based on these literature reports, one may select 
the value of flexural rigidity for a toroidal micelle. Us¬ 
ing such a range of representative values for the flexural 
rigidity, one can fit the theoretical size distribution $tot 
in (43) to the histogram of the experimental data to find 
the interval of values for tpo ■ 


nor intrinsic torsion. For a wormlike micelle consisting 
of more than one type of surfactant molecules, the pres¬ 
ence of molecules with different head-group or tail con¬ 
formations might lead to intrinsic curvature or intrinsic 
torsion. According to our modeling assumptions, such 
effects are ruled out. 

Using the derived free-energy density, the net free- 
energy of either closed or open wormlike micelles were 
derived. Whereas the total free-energy of a closed worm¬ 
like micelle is obtained simply by integrating its free- 
energy density over its centerline, the total free-energy of 
an open wormlike micelle includes two additional contri¬ 
butions from its end-caps. 

The derived model was applied on a special case of 
open wormlike micelles whose centerline has the shape of 
a circular arc. The model was also applied on a toruslike 
micelle. In both cases the minimum point corresponding 
to the equilibrium state was observed. 

The developed model was applied to toroidal worm¬ 
like micelles, and the results were compared with those 
available in the literature. The range of the minimum 
free-energy of a toroidal micelle was found to be on the 
same order of magnitude as the predictions previously 
reported. Further, our theory predicted that the equilib¬ 
rium state of toroidal micelles is stable for positive values 
of flexural rigidity, which is analogous to the prior predic¬ 
tions. Applying necessary statistical-thermodynamical 
concepts following previous studies, the theoretical size 
distribution of a toroidal micelle was found in terms of 
the derived free-energy. Consistent with the previous ob¬ 
servations, it was found that the size distribution of the 
toroidal micelle decreases by increasing the value of its 
flexural rigidity. Fitting the theoretical size distribution 
to the experimental data, the intervals of values for ma¬ 
terial parameters can be found. 

The derived model is capable of quantifying observa¬ 
tions of wormlike micelles. Comparison of the results 
with previous studies on wormlike micelles available in 
the literature confirms the accuracy of the model. 


VI. CONCLUSION 
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VIII. APPENDIX 


In this section, the details of the expansion of the 
right-hand side of (19) is carried out. On perform¬ 
ing the change of variables t — to = 6s and defining 
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r(s) = x(to) — x(io + < 5 s), ( 19 ) becomes 


^(L-tc)/S p2-ir p2-ir 

^ = si II 


l-to/S Jo Jo 


<52 


r(s) • e{to + Ss, 7]), d(to, 9 ) ■ e(to + Ss, r/)J 
f(to)f(to+Ss)d9d7)ds. ( 46 ) 


We wish to expend the right-hand side of ( 46 1 in powers 
of S neglecting terms of o(S^). The following abbrevia¬ 
tions are used: 


n:=n(to), t := t(to), b := b(to), 

f:=f(to), f':=f'(to), f":=f"(Q. 


( 47 ) 


To achieve the expansion of the right-hand side of ( 46 1, 
the following expansions are used: 

x(to + so) = x(to) + H-)n -I - —ktd 

zb 6 

-I- (s (5 - + o(S^), 

b 

2x2 2 x2 

n(to + sS) = (1 - 

2 x2 

— (sSk + + o(S^), 

2 x2 2 x2 

b(to + sS) = -(sSt + ^-^T')n -I- (1 - *-^r^)b 

s2A2 

-I-+ o(S^), 

/(to + sS) = f + sSf' -\ -+ o(( 5 ^). ( 48 ) 

By neglecting the higher-order terms and performing the 
dot products, the arguments of fl in the right-hand side 
of equation ( 46 1 become 


|x(to) - x(to + S5)|2 

P 


= s^ -b Aid'^s'^ + o{P), 


(x(to) - x(to + sS)) ■ d(to, 9 ) = A2S^s^ + o(S^), 
(x(to) - x(to + sS)) ■ e(to + sS, rj) = AsS'^s^ + o(( 5 ^), 
d(to, 9) ■ e{to + s6,r]) = A + A^Ss + A^d'^s^ + o((5^), 


( 49 ) 


where the coefficients A1-A5 in ( 49 ) are expressed as 
A = cos(0 —77), ~ ^2 = —^Kcosd, 

A3 =-K cos 77, A4 = T sin(0 — 77), 

A5 = - (r' sin (0 — T]) — cos 9 cos 77 — cos(d — 77)). 

( 50 ) 


The identities, t' = Kn, n' = —Kt + rb, and b' = —rn, 
called the Frenet-Serret formulas, which express the 
derivatives of t, n, and b with respect to the arclength 
s in terms of the Frenet frame {t,n,b} of C are used 
in ( 48 1.It should be noted that although the quantities 


Ai through A5 depend on 9 and 77, this dependence will 
not be denoted explicitly. Finally, put 


fl(s, 6 », 77 ) 

^,*(s,6»,77) 

(s, 9 ^ T]) 


:= ^(s^, 0, 0, cos(0 — 77)), 

:= il,j(s^,0,0,cos(6> - 77)), 

:= h,u (s^,0, 0, cos(6» - 77)), 

7 e { 1 , 2 , 3 , 4 }. ( 51 ) 


Using the expansions ( 49 ) and the abbreviations ( 51 ), 
while knowing <C and ^ and using ( 16 ), 


the right-hand side of ( 46 ) is, on neglecting terms pro¬ 
portional to and higher, 

/*27r /*27r 

= S I I I Cl{s,9,ri)f'^ d9dr]ds 


i-Uo Jo 


277 /‘27r 


1, 


^ / / 7 i^is,9,r])ff"s^ d 9 dr)ds 

J-iJo Jo ^ 

ct ,.27r,.27r 


- ^^i{s,9,r])s^ 


f i cos 9 cos 77114(5, 9, 77) 


d 9 


drj ds| 


2-17 p2-K f2 „2 


+ r 


— cos(0 — 77) 114(5,0,77) d6*d77d5| 


sin^(6> - 77)1144(5,0,77) 


( 52 ) 


which is the free-energy density in terms of the curvature 
and torsion. The final form of the free-energy density for 
a wormlike micelle is thus of the form 


where 

i’o = s 


4> = ljjo + kilP -I- fc 2 T^ 


p2TT p277 

I Ct{s, 9 ,ri)f'^ d 9 drjds 

-e Jo Jo 

•i n2TV r'iTT 


( 53 ) 


/t p^TT pZTZ -| 

J J -^^{3,9,7])/fs'^ d9dr]ds, 

j-t j' 2 -K j- 2 Tr _ r 2„2 . 

ki=S^ / / ———(6 cos 0 COS 77114(5, d, 77) 

J-eJo Jo V 


l-iJO Jo 
- 1 - 121 ( 5 , 0 , 77 ) 5 ^^ d0d77ds, 

f.27r p2-!r r2 „2 


fZnyZ17r2 2 

k2 = S^J J J ^-(^sin2(0- 77 ) 1144 ( 5 ,0,?7) 

— cos (0 — 77) 114(5,0,77)') d0d77d5. ( 54 ) 
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